Resonant scattering of ultracold atoms in low dimensions 
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Low energy scattering amplitudes for two atoms in one- and two-dimensional atomic wave guides 
are derived for short range isotropic and resonant interactions in high partial wave channels. Taking 
into account the finite width of the resonance which was neglected in previous works is shown to 
have important implications in the properties of the confinement induced resonances. For spin 
polarized fermions in quasi- ID wave guides it imposes a strong constraint on the atomic density for 
achieving the Fermi Tonks Girardeau gas. For a planar wave guide, the charateristics of the 2D 
induced scattering resonances in p- and d-wave are determined as a function of the 3D scattering 
parameters and of the wave guide frequency. 
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Recent experimental progress in degenerate atomic 
ases make possible accurate studies of quasi-one (ID) 
IBS II BM and quasi-two dimensionnal (2D) 0,0,11 
configurations. These systems are interesting in view of 
future applications involving coherent manipulation of 
matter waves and can be used also for studying generic 
phenomenona in low dimensions [H, HI]. One major inter- 
est of atomic gases is the precise knowledge and exper- 
imental control of the low energy interatomic collisions: 
the effective two-body interaction can be tuned with a 
Feshbach resonance |9| by applying an external magnetic 
field and/or with a Confined Induced Resonance (CIR) 
by varying the extension of the trap in the tight trans- 
verse direction [l(| HH, E3] ■ These techniques open pos- 
sibilities for achieving new types of strongly correlated 
quantum systems. For example, thanks to s-wave col- 
lisional properties in quasi-lD wave guides [Iol |. it has 
been possible to observe the so-called Tonks Girardeau 
(TG) gas where ID hard-core bosons can be mapped 
onto a system of non interacting fermions [H, H|. For 
quasi-lD spin polarized fermions an analogous regime: 
the Fermionic Tonks Girardeau (FTG) gas is subject to 
intensive studies [3, fl6| . In this case, as a con- 
sequence of the Pauli exclusion principle, atoms interact 
predominantly in the p-wave channel and it has been pre- 
dicted that the strongly interacting ID polarized Fermi 
gas can be mapped onto a non interacting Bose gas. 

In this letter, it is shown that the finite width of the 
high partial waves resonances is an essential feature in 
the properties of the corresponding CIR. As a conse- 
quence, the resonant fully polarized Fermi atomic gas in 
quasi-lD traps can reach the FTG regime in a very dilute 
limit only (which means that for realistic trap parame- 
ters, only few particles systems can undergo this regime) 
and is more generally described by a narrow BEC-BCS 
crossover. Violation of the fermion-boson mapping theo- 
rem [TtJ opens relevant issues on the quasi-lD resonant 
polarized Fermi gas such as the equation of state and the 
shift and damping of the collective modes which are only 
known at the FTG limit [l^|. Motivated by the exciting 
predictions of exotic superfluid phases [H, Qj| with non- 



Abelian statisitics, the p-wave scattering amplitude in 
quasi-2D atomic wave guide induced by a resonant inter- 
action is derived. The effective range parameter which 
were neglected in Ref.[2(| appears essential for the de- 
termination of the 2D low energy scattering parameters. 
Due to the possibility of achieving d-wave resonances [2l[ 
in 3D systems, the quasi-2D d-wave scattering amplitude 
is also derived and the characteristics of this new CIR 
are depicted. 

In the following, the true interatomic forces are sup- 
posed to be short range and isotropic. Moreover they are 
considered in the neighborhood of a resonance. The 3D 
low energy collisional properties of two atoms are then 
parameterized in each partial wave (I > 0) by the follow- 
ing phase shift Si (q) : 
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in Eq.([T]) q is the relative momentum of the two colliding 
atoms, and the resonant regime in the /-partial wave is 
achieved for \wi\ ^> R 2l+1 , where R is the characteristic 
radius of the pairwise potential (wq is the usual scatter- 
ing length and w± is the scattering volume) . The effective 
range parameter ai is linked to the width of the reso- 
nance: a large (small) value of aiR 21 ^ 1 corresponds to 
a narrow (relatively broad) resonance. In the neighbor- 
hood of the resonance and for I > 1, it has been shown 
recently that for a finite range potential [22|, [H, [24[ : 
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Thus ai\wi\ ^> R 2 , and —ti 2 / (2fj,aiwi) is a low energy 
scale (p is the reduced mass) which is for large and posi- 
tive values of wi , nothing but the shallow bound state en- 
ergy (and the quasi-bound state energy for large and neg- 
ative wi). Hence, unlike s-wave broad resonances where 
cvo = O(R) can be neglected in the low energy limit, for 
I > the effective range parameter ai which depends on 
the specific resonance considered, is an essential parame- 
ter involved in the low energy properties. In experiments 
the resonant regime can be achieved by using a Feshbach 
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resonance [9fl on the spin degree of freedom. An exter- 
nal magnetic field B modifies the detuning between an 
open and a closed channel of the two-body system, and in 
the vicinity of the resonance in the partial wave I where 
B ~ Bq, cx —(B — Bq) while ai is almost constant. 

Due to the short range character of interatomic forces, 
we use the so-called zero range approximation. Hence, 
the pairwise interaction between particles is replaced by a 
source term of vanishing range e — > in the Schrodinger's 
equation (it is assumed that the characteristic lengths of 
the trapping potential are ^> R). Moreover, the trapping 
potentials considered in this letter are always harmonic, 
so that the center of mass and relative coordinates are 
decoupled each from the other. For a binary collision 
in the center of mass frame at energy E = H 2 q 2 / Lb, the 
Schrodinger's equation can be transformed into an inte- 
gral equation for the wave function 
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where TLq is the free Hamiltonian which includes the 
external potential, and l^o) belongs to the kernel of 
Tio — E. In Eq.Q, the source term is introduced in 
the momentum representation: S e (k) = exp(— fe 2 e 2 /4) 
is the Fourier transform of a normalized Gaussian 
weight having a vanishing range e (25l |. IZi and Si^ are 
Symmetric Trace Free (STF) tensors of rank /, and the 
notation (IZrSi.u) means a contraction between the 
two tensors. The tensors Si^ are eigenfunctions of the 
momentum operator and they appear in a standard 
multipolar expansion in Cartesians coordinates [22T |: 
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(-l) l k l+1 (d km d k0 ...)k- 1 /{2l-l)\\, where 
{k a } are the Cartesian components of the vector k 
(a € {x, y, z}). Eq.© shows explicitly that the interact- 
ing wave function |*} is the superposition of a regular 
solution |*o) and of an irregular part generated by the 
source term (27| . The correct asymptotic behavior of 
the wave function (for relative coordinates |r| ^> R) is 
obtained by the determination of the tensors IZi which 
fixes the balance between the regular and irregular 
solutions. For this purpose, contact conditions (for 
r — ► 0) are imposed in each partial wave and are such 
that without external potential the phase shifts in a 
scattering process coincide exactly with the first two 
terms in the right hand side of Eq.|T|) (22J. In the 
momentum representation, the contact conditions can 
be written as: 
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where af 1 (q) — (itf; + mq 2 ). In Eq.Q, Reg e _ +0 means 
the regular part of the integral obtained when the formal 
range e is set to zero. This way, the source term in Eq.Q 
is itself a functional of |*), and a closed equation for IZi is 
obtained by combination of Eqs. (l3l4p . In the next parts, 
the scattering problem for p- or ci- wave channels is solved 
successively in ID and 2D harmonic wave guides. 



Linear atomic wave guide. — Two atoms are confined 
in a two-dimensional harmonic trap while they move 
freely along the third direction (z). In the center of mass 
frame, the non interacting Hamiltonian is: 



h 2 

Ho = -^A r 
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where r = ze z + p are the relative coordinates. For a 
scattering process at energy E = h 2 q 2 /2/x, the state |^o) 
in Eq.© is: (k|* > = (2n)S(k z ~ q) (k 2D |0 O o), where 
\4>oo) is the ground state of the 2D harmonic oscillator 
in Eq.©, and k = k 2 n + k z e z . Hereafter, the system is 
in the monomode regime (E < 2frui±). At large distances 
\z\ ^> a± (where a± = ^/h/nu>± ^> R), the wave function 
factorizes as (r|*) ~ {p\4'oo}' 1 Pid( z ), and for E > 0: 



*1>id(z) = exp(igz) + exp(zg|z|) 



sign(z)/° 



.(6) 



In this situation, the system which is frozen along the 
transverse direction, can be considered as quasi-lD. The 
even (odd) scattering amplitude / even (/ odd ) results from 
the asymptotic contributions (z a±) in the different 
even (odd) 3D partial waves. Quasi-lD scattering in the 
s-wave channel has been thoroughly studied [10( and in 
the following we consider quasi-lD scattering of two spin 
polarized fermions in the vicinity of a p-wave resonance. 
Thus, one can restrict the source term in Eq.Q to the 
I = 1 contribution, where <Si,k = k/fc, and the choice |*o) 
imposes that IZi = Pid&z- Computation of P\d from 
Eq. (|4]) is much simpler in the domain of negative energy 
E < and using Eq.© one can show that: 
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where r = Ej2hw^ < and P.f. denotes the "partie 
finie de Hadamard" of the integral (28[. Interest- 
ingly, one recognizes the Riemann-Hurwitz Zeta func- 
tion £h-(— 1/2, — r) in the regularized integral of Eq.([7]). 
Hence, by analytic continuation in the domain r > 0, one 
obtains a simple expression of the scattering amplitude, 
and in the low energy limit qa± <C 1: 



f° dd = 2bTrP 1D <f>* (0)~-iq[ - +iq + q 2 i p ) , (8) 
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with the odd-wave scattering length l p and the effective 
range £ p given respectively by: 



6aj_ 



a>ia 2 , 

and £ p = ^±. (9) 



Expression of l p in Eq.© coincides with the result in 
Ref. [l2j , while the effective range £ p and its physical im- 
plications was (up to our knowledge) not studied in pre- 
vious works. However, £ p is in general not negligible for 
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all gCa^ 1 and crucially depends on the width of the 
3D resonance and on a±. Indeed, the inequality £ p 3> a± 
is likely to occur as a consequence of Eq.((2|) (a±R > 1) 
and also from the condition a± S> R which is needed 
from the hypothesis that scattering processes are studied 
at collisional energies <§; h 2 /fj,R 2 . The regime £ p ~ a± 
can be reached only for an extreme transverse confine- 
ment (a± ~ lOi?) and for the broadest p-wave resonances 
with a±R r-j 1. In actual experiments, two p-wave reso- 
nances are used: the resonance for 6 Li atoms (R ~ 3 nm) 
at B ~ 215 G with a\R ~ 5 [2!| and the one at 
Bq ~ 198.8 G for 40 K atoms (i? ~ 7 nm) where a x R ~ 3 
[301 ] . For Z p > the scattering amplitude in Eq.® has a 
pole at q = in with k = (—1 + yl + A^ p /l p )/2^ p > giv- 
ing a shallow bound state energy at (— h 2 n 2 /2ji). For 



large and positive values of L 3> £ 



Z^ 1 and the 

bound state has a vanishing energy (i.e. <C hw±) [3~j| . 
The scattering amplitude in Eq. (jSj) can be obtained from 
a ID effective theory where the wave function ipiD in 
Eq.© solves the non interacting Schrodinger equation 
and satisfies the following contact condition: 

Hm (J- + d z Cpd^j ^(z) = 0, (10) 

where tpi^(z) — [ipn)(z) — ipii)(—z)]/2 is the projection 
of the wave function onto its odd component [32J. This 
approach can be generalized for few- and many-body sys- 
tems by imposing the contact condition in Eq. (|10p for 
each pair of interacting particles. Without performing 
these calculations which are clearly beyond the scope 
of this letter, it is of importance to determine the con- 
ditions such that the fully polarized fermionic gas can 
reach the FTG regime or equivalently can be mapped 
onto non interacting ID bosons [T3|. Eq.JTU]) implies that 
this mapping is possible only at resonance (\l p \ = oo) 
and also if the effective range is negligible. The reso- 
nant condition is easily obtained by using a Feshbach 
resonance. However, the momentum distribution of the 
FTG gas has a large tail given by a Lorentzian of width 
An, where n is the ID atomic density [HI], thus £ p can 
be neglected only in the dilute limit: n£ p -C 1. For N 
atoms trapped in a strongly anisotropic trap with a weak 
harmonic confinement along the z-direction (atomic fre- 
quency lo z <C u>± and axial length a z = y/H/ uoj z ^> a±), 
in the FTG regime one has n ~ N / a z and the condi- 
tion n£p -C 1 gives N -C 6a z /(aia 2 j_). For example, con- 
sidering 40 K atoms in a highly anisotropic trap with 
u>± = 2ir x 70 kHz and uj z = 2ir x 10 Hz, the FTG regime 
is obtained for TV -C 14, which makes sense only for few- 
body configurations (for 6 Li atoms with the same trap 
parameters, £ p is of the order of a z and the mapping to 
a non interacting Bose system is a poor approximation 
even for the two-body ground state). Interestingly, if the 
condition n£ p < 1 is not satisfied, the scale invariance 
in the linear wave guide is broken at low energy, hence 
the time-dependent many-body ansatz in Ref.[lqj is no 
more an exact solution. To conclude this part, excepted 
specific configurations, the criterium n£ p <C 1 is in gen- 



eral not verified and the density corrections to the FTG 
properties are important issues for understanding the res- 
onant gas. Moreover, following the reasoning of the box 
model in Refs.[23|, Hj], one expects that by varying l p 
from large and positive values to large and negative val- 
ues, the system experiences a narrow p-wave BEC-BCS 
cross-over, where the composite bosons in the dilute BEC 
phase (lp > 0) are dimers populating the quasi-lD shal- 
low two-body bound state. 

Planar atomic wave guide. — High partial wave super- 
fluidity in quasi-2D geometries is interesting for its links 
with condensed matter physics like for example high-T c 
superconductivity and the possible applications for quan- 
tum computing [1|| EH ■ These studies motivate a close 
investigation of quasi-2D scattering properties. In this 
geometry, the two colliding atoms are confined in a pla- 
nar harmonic trap along the z-direction and move freely 
in the (xy) plane. The free Hamiltonian reads: 



H = -— A r 
2u 
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The homogeneous solution corresponding to a scat- 
tering process at energy E — h 2 q 2 /2u is: (k|^o) = 
(27r) 2 ^(k 2 £) — q) (k z \4> ), where \<f>o) is the ground state of 
the ID harmonic oscillator in Eq. (fTTj) . In the monomode 
regime (E < hu z ), the system can be considered as quasi- 
2D, and for large interparticle separation (p » a z ), the 
wave function factorizes: \&(r) = </> (z)tp2D (p)- The 2D 
partial scattering amplitudes /I™' = /[~ m l can be defined 
by the following expansion: 



MP) = e<*»-"> 
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f^H^(qp)eJ" ne \ (12) 



where 9 — n/2 + Z(p, q) and Hm is the Hankel's 
function. Using Eqs. I|3I4I12[) the scattering ampli- 
tude for two atoms interacting in the s-wave chan- 
nel is: /M = Ai:[a z ^ / a Q (q) + J (r + z0+)]-\ with 
r = E/2hw z , J (r) = Iu(-B/(2ttt)) + En=x^(n/(n - 
r))(2n - l)!!/(2n)!! H, and B ~ 0.9049 [36|. In the 
p-wave channel, Eq.® implies the existence of a small 
parameter r/i — (aia 2 ) _1 <C 1 which plays a central role 
in the scattering properties [13] • For q = qe x , the p-wave 
regular part in Eq. (J3J) is related to the 2D p-wave scatter- 
ing amplitude /W by K x = -/W(«)gx/[27r#S(0)]. Using 
the contact condition in Eq.([5]) where S\^. — k/k gives: 



/W( g )-6Vl0o(O)| 5 



1 



6Ji(r 



i0+) 



n -1 



ai(q) 



(13) 



where |0o(O)| 2 = l/(a z y/ir). In the low energy limit 
|r| -> then J x (t) = J x (0) + r ln(-eS/2vrr) + 0(r 2 ), 
where J x (0) ~ -5.4722 x 10~ 2 35], and the logarith- 
mic term ensures the unitarity condition in Eq. (|13p . 
A similar calculation can be done for atoms interact- 
ing resonantly in the 3D <i-wave channel. In this case, 
the small parameter related to the resonance width is 
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T]2 = a 2 a z 3 <C 1. For q = qe x , the I = 2 regular part is 
given by: k 2 (TZ 2 -S 2ik ) = S 2D (2k 2 z - k 2 D ) + D 2D (k 2 x - fc 2 ,), 
where S2.D contributes to the m — 2D partial wave 
channel and D 2 d contributes to the \m\ = 2 channel (2D- 
(i-wave). The scattering amplitude in the \m\ = 2 channel 
is /[21(g) = -2Trq 2 (j)* (0)D 2D , and using Eq.®: 



/ [21 (9) 



157r g 4 |0 o (O)| ; 



1 



60 J 2 (r + i0+) 



(14) 

where J 2 (r) = J 2 (0) + rJ^O) + 6>(r 2 ) f or r - » 0, and 
J 2 (0) ~ -2.2752 x 10~ 2 35]. From Eqs.fHIl]), several 
conclusions can be drawn on the quasi-2D scattering in 
m = 1 and m — 2 partial wave channels: i) At small 
collisional energy (q — ► 0), f^ m \q) oc q 2m /(w" 1 — w^ 1 ) 



where 



5.39 x a?, and 



~ 1.3 x a%. Hence, 
the 2D confinement induces a shift w* n in the general- 
ized scattering length w rn , which grows as the confine- 
ment increases, ii) By modifying the trap frequency 
and/or the external magnetic field it is possible to drive 
the system from a regime with a 2D m-wave quasi- 
bound state (w^ 1 < w^ 1 ) to a regime with a shallow 



bound state (0 < w m < w^), in both cases the (quasi- 
) bound state energy is Ef, ~ —h 2 (w^ — w*^ 1 ) / '(2/xa m ) 
|38| . thus having an expression similar to the 3D case: 
Efr D ~ —h 2 / (2^iwiai). iii) Resonance in the scattering 
cross section a — |/ m | 2 /4g occurs at a collisional energy 
E Et, that is only for positive values of E^: in pres- 
ence of a quasi-bound state. The resonance width is given 
by AE/E b oc ^(Eb/fkVz)" 1 ^ 1 <C 1. iv) Consequently, a t 
appears as a crucial parameter for the low energy scatter- 
ing properties. For example, by neglecting a\ (rji = oo) 
in Ref . [201 ] , the resonant collisional energy and resonance 
width of the p-wave CIR are found with other order 
of magnitudes, v) Following the reasoning in Ref. [13], 
these scattering properties open the possibility of observ- 
ing high partial waves BEC-BCS transitions in quasi-2D 
fermionic gases by varying (if" 1 — w^ 1 ) from positive 
to negative values. 
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